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Abstract 

We present a method which enables exact solutions to be found for fiat 
homogeneous and isotropic scalar-tensor cosmologies with an arbitrary cj(<1?) 
function and satisfying the general perfect fluid state equation P = (7 — l)pc 2 . 
This method has been used to analyze a wide range of asymptotic analytical 
solutions at early and late times for different epochs in the cosmic history: 
false vacuum inflationary models, vacuum and radiation-dominated models, 
and matter-dominated models. We also describe the qualitative behavior of 
models at intermediary times and give exact solutions at any time for some 
particular scalar-tensor theories 

PACS number(s): 04.50.+h, 04.80,+z, 98.80.Cq, 98.80.Hw 
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I. INTRODUCTION 



The simplest generalization of Einstein's theory of gravity are scalar-tensor theories, 
which contain the metric tensor g^, and an additional scalar field <fi. The strength of the 
coupling between gravity and the scalar field is determined by an arbitrary coupling function 
u)(4>). Although these theories have a long history M, they have received a renewed interest 
both in cosmology and particle physics mainly because of three reasons. First, scalar-tensor 
theories are important in inflationary cosmology since they provide a way of exiting the 
inflationary epoch in a non- fine-tuned way ||. Second, most of the recent attempts at 
unified models of fundamental interactions predict the existence of scalar partners to the 
tensor gravity of General Relativity (GR). This is the case of various Kaluza-Klein theories 
@, as well as supersymmetric theories with extra dimensions and superstring theories [|J, 
where scalar-tensor theories appear naturally as a low-energy limit. Finally, these theories 
can satisfy all the weak-field solar-system experiments and other present observations || 
to arbitrary accuracy, but they still diverge from GR in the strong limit. Thus, scalar- 
tensor theories provide an important framework for comparison with results of gravitational 
experiments which could take place in the near future (e.g., LIGOS or VIRGO 0). 

Analytical or numerical solutions for the scalar-tensor cosmological models are well 
known in the framework of some particular theories proposed in the literature. This is 
the case of Brans- Dicke's theory ||, Barker's constant-G theory ||, Bekenstein's variable 
rest mass theory []10 |, or Schmidt-Greiner-Heinz-Muller's theory fll |. In the last years, a 
considerable effort has been devoted to investigate the cosmological models in more general 
scalar-tensor theories. For example, Burd and Coley [|12|] have used a dynamical system 
treatment to analyze the qualitative behavior of models where a constant (Brans-Dicke) 
coupling function is perturbed by a slight dependence on the scalar field. 

A great improvement in the search of solutions of the scalar-tensor field equations has 
recently arisen in the form of methods which allow analytical integration through suitable 
changes of variables. This approach was first presented by Barrow who generalized the 
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method by Lorentz-Petzold |14| to solve the vacuum and radiation-dominated cosmological 
equations of any scalar-tensor theory By using and extending this method, Serna & Alimi 
fl"5| have performed an exhaustive study of the early-time behavior during the radiation- 
dominated epoch of scalar-tensor theories where u(4>) is a monotonic, but arbitrary, function 
of 0. 

A further step in the search of exact scalar-tensor cosmological solutions was provided by 



Barrow, Mimoso and others who developed a method to solve the Friedmann-Robertson- 
Walker field equations in models with a perfect fluid satisfying the equation of state P = 
(7 — l)pc 2 (with 7 a constant and < 7 < 2). In this method, solutions are found through a 
generating function g which is indirectly related to the coupling function u;(0). A systematic 
study of the asymptotic behavior of scalar-tensor theories, although possible, is then difficult 
to be performed since the connection between g and uj(4>) is only known after completely 
solving the field equations. Exact solutions during all the epochs of cosmological interest 



have been nevertheless derived for some classes of scalar-tensor theories p6| . |T7 [. 

We will present in this paper an alternative method, also based on a suitable change of 
variables, which enables exact solutions to be found for homogeneous and isotropic scalar- 
tensor cosmologies with an arbitrary uj(<p) function. Unlike the previous methods (except for 
some studies of the late time behavior of scalar-tensor theories [18]] ), the generating function 
will be now the time dependence of the coupling function itself, what will allow for an easier 
analysis of the early and late time asymptotic solutions during all the main epochs in the 
universe evolution. 

The paper is arranged as follows. We begin outlining the scalar-tensor theories (Sect. |1D 
and we then present a method to build up homogeneous and isotropic cosmological models 
in their framework (Sect. PTT|) . This method is then applied to analyze the early (Sects. [TV] 
and |V|) and late time solutions (Sect. [VT] ) for a wide range of theories. Finally, conclusions 
and a summary of our results are given in Sec. |V11. 
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II. SCALAR-TENSOR COSMOLOGIES 



A. Field equations 



The most general action describing a massless scalar-tensor theory of gravitation is [IS 



s = ~Ll m ~ ^f-M")Vgd^ + s M , (i) 

where TZ is the curvature scalar of the metrics g^, g = det(g fJiU ), and is a scalar field. 
The strength of the coupling between gravity and the scalar field is determined by the 
arbitrary function (J ((f)), usually termed the coupling function. Each specific choice of the 
u)((f>) function defines a particular scalar-tensor theory The simplest case is Brans-Dicke's 
theory, where u is a constant. 

The variation of Eq. ([[]) with respect to g^ v and (f> leads to the field equations: 



1 „ 8vr m u , , , 1 
T 



ft^u - 7^9^ = — — - —((t)^(t),u - -9llv4>,k4>' K ) 



~((f>,w- 9^4>) (2a) 
(3 + 2uj) \3<f> = 8nT - c/0, K 0' K (2b) 

which satisfy the usual conservation law 

Tp = 0, (3) 

where T^ u is the energy-momentum tensor, uj' denotes duj/d(f> and □</> = g^ u (f) tfl - u . 

B. Cosmological models 

In order to build up cosmological models in the framework of scalar-tensor gravity the- 
ories, we consider a homogeneous and isotropic universe. The line-element has then a 
Robert son- Walker form: 



ds 2 = -c 2 dt 2 + a 2 (t) 



(l ' fl + r 2 dQ 2 



l-Kr 



(4) 



4 



and the energy-momentum tensor corresponds to that of a perfect fluid 



= (p + P/c 2 )u IM u u + Pg i 



(5) 



where K = 0, ±1, a(t) is the scale factor, p and P are the energy-mass density and pressure, 
respectively, and is the 4-velocity of the fluid. 
By writing the equation of state as 



P = ( 7 - l)pc 2 



(0 < 7 < 2) 



and assuming a flat Universe (K = 0), the field equations (Q) become 



8nGp 
3$ 



2$ 



(3 + 2u) 



a $ 
a + 2$ 



CO 



a 

a 3 + 2uj 



87rGp(4 - 37) 



(6) 

(7a) 
(7b) 



(3 + 2lu) 

where G is Newton's gravitational constant, $ = G<f>, and dots mean time derivatives. In 
addition, the conservation equation @ implies 

-37 



p oc a 



(8) 



III. INTEGRATION METHOD 

The main difficulty for finding complete or asymptotic solutions of the scalar-tensor 
cosmological equations is the presence of a source term in the right hand side of the wave 
equation ([7b]) . The only case in which these equations can be solved for universes of all 
curvatures is 7 = 4/3 (vacuum and radiation-dominated models), where Eq. ([TG|) becomes 
sourceless [0. For general perfect fluids (7 7^ 4/3), Eqs. (0) can be solved only for flat 



universes and using 'indirect' methods [0 where a particular scalar-tensor theory is defined 
through some generating function, g, instead of the coupling function cj($). 

We will now present a 'semi-indirect' method where the coupling function oj is used as 
the generating function defining each particular theory. However, unlike the 'direct' method 
available for radiation-dominated (7 = 4/3) epochs, the form of to>($) as a function of $ is 
not initially known, but only its dependence on a 'time' variable r. 



A. Integration method 



In order to solve Eqs. (H), we introduce a new time variable r 

dr = VTa^-^dt 



(9) 



and two dynamical variables 



y = a 3 ( 2 -^$y I 3 + 2u | /3 
x = a 3(2 ~ 7) $ 



(10) 
(11) 



where primes denote differentiation with respect to r and V = 87rGpa 37 /3 = constant. 
Using these variables, the field equations (^) reduce to 



y = (4 - 3 7 ) 



a 



3 + 2w | /3 



(12) 



and 



x + cn/ 2 /4 

where cr = sign(3 + 2ui). 

By introducing the function 



(6 - 37) 2 



/ , (y 

x + a — 



2V 



(13) 



f = x + (xy 2 ^, 



Eq. (jl3|) can be written in a more compact way as 



(14) 



(/ 2 )' = 3(2- 7 )/ 

Integration of Eqs. (|12D and (|1^) is then straightforward. The solutions are 

dr 



(15) 



y = o-(4 - 37) 



/=-(2- 7 )t + /o 



3 + 2cj I /3 



+ 2/o 



(16) 
(17) 



where and /o are integration constants. We see from Eqs. ([II]) and ( |T7D that the particular 
cases 4 — 37 = and 3(2 — 7) = imply y =constant and / =constant, respectively. In the 
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first case, complete exact solutions can be easily found from Eq. ( |T0"D because the second 
integration implied by Eq. ( |T6D is not needed. Since the 3 + 2uj function, as well as the 
constant factors (4 — 37) and 3(2 — 7) will appear several times throughout this paper, we 
will denote 

W(t) = (3 + 2w)/3 (18) 

71 = 3(2 - 7) (19) 

72 = 4 - 3 7 (20) 

Constants 71 and 72 are known at any epoch in the universe evolution. A hypothetical 
false- vacuum era would correspond to 7 = (71 = 6, 72 = 4), while the vacuum (r = 0) 
and radiation- dominated (r > 0) eras correspond to 7 = 4/3 (71 = 2, 72 = 0). Finally, the 
matter-dominated era is defined by 7 = 1 (71 = 3, 72 = 1). 

We will see now how Eqs. (|9"D-(|T7D enable exact solutions to be found for the scalar-tensor 
field equations (|7|). To that end, we will consider separately the cases 7 < 4/3 (72 > 0) and 
7 = 4/3 (72 = 0). 

1. Case 7 < 4/3 (general perfect fluids) 
When 72 > 0, we can use Eqs. (|i0l)-(|i~2j) to write 



$ 72 X 



(21) 

By integrating this equation, and using Eq. (|I~4"D, we find 

$ = $ 1 (e J / |x|) 2/72 ( 22) 

where sign($) = sign(x) (see Eq. [11]), and 

J = 71 [ T -dr (23) 



Tl X 



with $1 and n being integration constants. 



From Eqs. ([TTJ) and (0), the scale factor is then given by 

a = ai (\x\ / e 2J/7l ) 1/72 (24) 

where a\ > is another integration constant. 

The a(t), $(£) and functions can be then obtained in the following way: i) we start 
by choosing a particular scalar-tensor theory, that is, by specifying 3 + 2uo as a function of 
r; m) we compute y(r) from Eq. (|I6|); in) x{t) is then obtained from Eqs. (fUD and (|i~7|); 
to) $(t) from Eq. (0); u) cj($) by inverting $(r) and using 3 + 2u(t); vi) a(r) from Eq. 
(pjj); and mi) r(t) from Eq. (|9|). Obviously, we assume that the corresponding expressions 
are all integrable and invertible. 

Once a(r) is known, we can find other functions of cosmo logical interest. For example, 
from Eq. (0) , the Hubble parameter H = a/a is given by 

H = Vfa 3 ^ ( -) . (25) 



Differentiating Eqs. QTTj ) and (H), and using Eqs. ( JTOj ) and (|T^) to eliminate $ and / 



derivatives, h = a'/a can be written in terms of /(t), y(r), and w(r) as 

(26) 



s [" 2^| 3 + 2w| /3_ 
Another function of interest is the speed-up factor £ = H/Hgr, which compares the 
expansion rate, H, to that obtained in the framework of GR at the same temperature, Hqr- 
To obtain this quantity, we use Eq. (p5|) and Hq R = 8nGp/3 = Ta" 37 , so that 

£ = a^' 2 a' = a^ /2 h (27) 



2. Case 7 = 4/3 (vacuum and radiation- dominated models) 



When 72 = 0, Eqs. (p6|)-(p7D hold but now with a constant y value (see Eq. |T2"D : 



y = yo 



(28) 
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Eqs. ( |TOD and fll4]) then become 

&a 2 = V3y \3 + 2uj\- 1/2 (29) 



(r + h) 2 - cry 2 / A (30) 



which coincide with Eqs. (15)-(16) of Serna & Alimi Jl5] (see also ref. [0). 

Using Eq. (|TTD to eliminate the scale factor in Eq. (|29"D, and integrating Eqs. (|29|)- (|30"|). 
we obtain (for r ^ 0): 

/|3 + 2u;| 1/2 f= (31) 

^{ttM) (*>°) 

2V3 [atan (^) - sign(A)7r/2j (a < 0) 



where A = yo/2. 

From a similar procedure (but using adr = c?t), it is straightforward to show that, when 
T = and a > 

J |3 + 2^| 1 / 2 ^ = v^sign(A)ln[|A|(T + /o)] (32) 

while no solutions exist for a < 0. 

Since a detailed study of the radiation-dominated asymptotic solutions has been per- 
formed in a previous paper |L5|], we will give here just a tabulated summary of results (see 



Table 1). 

B. Setting of constants 

Since the origin of times can be chosen according to any arbitrary criterion, an analysis 
of asymptotic solutions at early times (r — > 0) would be meaningless if we do not clearly 
specify what means r = for us. In this paper, we will set the origin of times by using 



a criterion similar to that of Serna & Alimi |Tj|. Using the notation of the present paper 



and assuming that the sign of 3 + 2uj never changes during the Universe evolution, such a 
criterion can be expressed as follows: 
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1. Models with 3 + 2uj > (a = +1) 



As we will see now, these models are characterized by the fact that x vanishes at least 
at one finite value of r. Consequently, we take the origin of times so that the most recent 



zero in x occurs at r = (x(0) = 0, and x > for any r > 0). We see from Eq. (|TT|) that 
such a criterion ensures two important features of our solutions: 1) a singular point either 
in the scale factor a or in the scalar field $ defines the origin of times, 2) both a and $ have 
positive values (strictly) for any r > 0. 

From Eq. flT4]), we see that condition x(0) = can be written as 

fo = <ry 2 o/* ■ (33) 

Throughout this paper, the subindex denotes value at r = 0. For a = +1 and any 
value of yo (positive, negative or vanishing), we then find that the fo value can always be 
chosen so that x(0) = 0. 

The condition x > for any r > implies that f > 0. As a matter of fact, if fo were 
negative, the function / = (7i/2)r + /o would vanish at a positive time r* = —2/0/71. Thus, 
x = f 2 — y 2 /4 = — y 2 /4 < at T* > 0, in contradiction with the requirement x > for any 
t > 0. By combining fo > with the condition (|33|), we then find that integration constants 
must be chosen as 

/o = M/2 (<r = +l). (34) 

Another important consequence of the requirement x > for any r > is that / > 
I y I o/2. As a matter of fact, when this requirement is combined with Eq. (0), one obtains 
that /(r) > \y{r)\/2. Since f = |z/o|/2, the last inequality can only be satisfied if / > |y|' /2. 
Consequently, from Eqs. ([12]) and (|17D, we find that 



^Wt. (if */»><>) (j = +i) (35) 

(if !/„ < 0) 



We note however that Eqs. (|33J) and (^) are just necessary, but not sufficient, conditions 
to satisfy our criterion for the origin of times. 
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2. Models with 3 + 2uj < (a = -I) 
When 3 + 2u < (i.e., o = — 1), the only choice satisfying Eq. (|3~3"D is fo = yo = 0. That 

is. 

/o = (if y = and a = -1) . (36) 

For a = —1 and non- vanishing (positive or negative) yo values, the condition ([^J) is not 
satisfied and x then never vanishes. In that case, we will see now that the scale factor has 
a non-zero minimum value and models are all non-singular. Consequently, we will take the 
origin of times so that the last stationary point of a occurs at r = 0. 

The existence of maxima or minima in the time evolution of the scale factor has been 



previously analyzed only for vacuum and radiation-dominated scalar-tensor cosmologies |L7 . 
The integration method presented in this paper also allows for an easy study of the existence 
of stationary points at any epoch of the universe evolution. 

The condition for the scale factor to contain a stationary point, a 1 = with a > 0, is 
equivalent to ho = 0, which leads to 

fo = - J !^= (if Vo + and a = -1) (37) 
2a/ W 



as the appropriate condition to take the origin of times so that a has initially a stationary 
point. 

We must note that the integration constants fo and yo must have finite values. Otherwise, 
solutions would diverge at any time. Consequently, Eq. ( pTf ) implies that the possibility 
Wq — > is excluded when y ^ and 3 + Iuj < 0. 



C. The coupling function 

1. Representation for early and late times 
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We will use the same representation for the coupling function as in our previous work 
, 20fl concerning radiation-dominated cosmo logical models. That is, 
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|3 + 2cu|= (3/A 2 )(A;+ (38) 

where 1/2 < 5 < 2 (to allow for convergence towards GR), A > 0, and k are arbitrary 
constants. This form gives an exact description for most of the particular scalar-tensor 
theories proposed in the literature and, in addition, it contains all the possible early behaviors 
of any theory where u;($) is a monotonic, but arbitrary, function of $ (see Appendix A). 

Barrow & Parsons |L7j have already analyzed a particular case {k = 0) of the coupling 
function given in Eq. (|3^). The k ^ term included here is nevertheless very interesting 
because it allows for models with $ — > oo and 3 + 2uj — > 3k /X 2 > as t — > 0. As shown 
by Serna & Alimi |p0| |, a wide class of such models are able to significantly deviate from 



GR during primordial nucleosynthesis and still predict the observed primordial abundance 
of light elements. In addition, such models allow for baryon densities much larger than in 
GR. 

We will just analyse the models where the asymptotic functional form of the scalar field 
$ is a power law of r: <J> oc r 13 or $ oc r~^, where (3 > 0. However as we will see later, this 
r-dependance for the $ field is found for most of the solutions. Consequently: 

(i) When $ = cr _/3 — > oo at r — > 0, we can approximate $ — 1 ~ $ = cr^ 13 (c =const.) 
and Eq. ( |3"8"D becomes 

|3 + 2cj| /3 = \W(t)\ = a 2 + br e (39) 

where b = (A 2 c a ) _1 , e = (35, and a 2 = k/X 2 denotes the | 3 + 2oj \ /3 value at r = 0. 

(ii) When $ = ct 13 — >• at r — > 0, we can approximate 1/(1 — $) 5 ~ 1 + 5$ = 1 + c^r^ 
and the coupling function is again given by Eq. ( |39|) but now with b = c5/X 2 , e = (3, and 
a 2 = (k+l)/X 2 . 



2. Intermediary times 



Although co>($) can be any arbitrary function of $, in order to discuss the behavior of 
scalar-tensor cosmologies at intermediary times we must impose some conditions on the 
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coupling function. 

Compatibility with solar-system experiments is only possible when the coupling function 
converges to the GR value (u — > oo) at late times. The simplest way of ensuring this 
condition is to impose that u(<&) is a monotonic (but arbitrary) function of $ so that it 
becomes infinity when $ — > 1. In addition, we chose the coupling function so that it implies 
a scalar field which is a monotonic function of r (i.e., the variable y of Eq. [16] does not 
become zero at a non- vanishing finite time). Consequently, the function W(t) will be also 
a monotonic function of time. 

An important consequence of this choice for the coupling function is that Eqs. (|34[) and 
(|35|) are now necessary and sufficient conditions to ensure the requirement x{r) > for any 
t > 0. As a matter of fact, if o = +1 and yo < 0, function y(r) must be negative at any 
time and, from Eq. fll2|), its absolute value \y(r)\ will be a monotonic decreasing function 
of t. Since Eq. (|T7|) implies that f(r) is always a monotonic increasing function, we deduce 
that x(t) = f 2 — y 2 /4: is also a monotonic increasing function of r. Consequently, x(0) = 
implies that x{r) > for any r > 0. 

On the other hand, if a = +1 and y$ > 0, the monotonic behavior of W(t), together with 



Eqs. (|g), QXTP and imply y'{r)/2 = l2 /(2^\W(r)\ < 7l /2 = f'(r), that is, x'{r) > 0. 



Consequently, x(0) = implies that x(r) > for any r > 0. 

IV. MODELS WITH 3 + 2uj > 

The analysis of models with 3 + 2uj > is specially simple because, in that case, the 
variable x can be written as (see Eq. [14]) 



x 



u(t)v(t) 



(40) 



with 



u = f + y/2 



1 i rz i \ i 



(41) 



v = f-y/2 




(42) 
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Here u = fo + y /2, v = f - y /2, while 



dr 



| 3 + 2lu I /3 

In terms of u and v, Eqs. (|22|)-(|2~4|) can be expressed as 



(43) 



n 

2 



a y 2 



2Ju 



2Jy 



(44) 
(45) 



with 



2 Jti u ' 2 J T1 v 



(46) 



At early times, the coupling function can be represented by Eq. (^) which, introduced 
into Eq. ([i"3|) , implies 



9(t) = 

and, hence, it and v are given by: 
If a > 0: 



-T 



1 - 



2a 2 (l+e) 



+ 



2t( 2 - £ )/ 2 
\/6(2-e) 



(« > 0) 

(a = 0) 



If a = 0: 



(47) 



(48) 
(49) 



, 7l , 72 (2-e)/2 

2 ^(2 - e) 

, 71 72 (2-e)/2 

v = v H r = r l ;/ 

2 V6(2 - e) 



(50) 
(51) 



The early-time evolution of £, a, and $ can be then obtained by using Eqs. (p7|), (f44])-([46"|), 
and approximating the it and v functions by the first non-vanishing term in Eqs. (|4"8])- 
(pT|). Furthermore, the integration constants appearing in this procedure must be chosen 



according to the criteria specified in section [I! IB 
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We will now analyze both the general features and early-time solutions of scalar-tensor 
models by considering separately the different possible values of constants a and dq. Both 
aspects will be illustrated by an example solved analytically by using Eqs. (|40|)-(|46|) or, in 
some cases, numerically 



A. Models with and y < ($'(0) < 0) 

1. General features at arbitrary times 



Since yo < 0, conditions established in Section |III C| imply that $(r) is a monotonic de- 
creasing function (y(r) < 0) and, therefore, the scale factor a(r) has a monotonic expansion 
without stationary points (see Eq. p6|) . In addition, the requirement x(0) — implies the 
existence of an initial singularity (a(0) = 0). 

These behaviors are shown in Figures 1, which present numerical solutions obtained by 
considering those theories defined, at any time, by W(r) — (p + qr) 3 (with a 2 = p 3 > 0, 
b = 3p 2 q > 0). 



2. Early-time solutions 



As shown in Sect. jTTBl, when y < 0, a necessary condition to have the latest singularity 



at t = is that integration constants satisfy fo = —yo/2 (see Eq. |34|). Then, as a first-order 
approximation, Eqs. fl48D and ( flDf) reduce to 



v = -yo 



(52) 
(53) 



Using Eqs. 



-m 



and (|9|), the early time behavior of models is then given by 



-2 -2 

$ OC T Jia+ "' 2 OC t 3a + 1 



(54a) 
(54b) 
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(2-37) + 71Q , , 3 7 a + 1 

£ oc r 2 (T2+7i") oc t + 2 3"+! (54c) 

3a + l 

£ oc r^i Q +^2 (54d) 

As expected, all models are singular (a(0) = 0) and $(0) — > +oo. Thus, we can consider 
$ ~| $ - 1 | at early-times so that, by identifying a = \/k/\ and e = 28/^xoc + 72), the 
coupling function (|3D| ) reduces to the form of Eq. ( ft£| ) . 

From Eq. ( |54cD , we also find that the early behavior of the speed-up factor depends on 
the a parameter: £ initially becomes zero (£(0) — > 0), finite (£(0) — > £0 > 0), or infinity 
(£(0) — ► +00) depending on whether a is smaller, equal or larger, respectively, than a critical 
value a c = (37 — 2)/^%. These three kinds of models are shown in Figures 1, where we also 
see that the case £(0) = has a non-monotonic behavior. 

B. Models with a ^ and y > ($'(0) > 0) 
1. General features at arbitrary times 



According to Eqs. (33) and (^), integration constants in these models must satisfy 
fo — 2/o/2. When this condition is introduced into Eq. (p6|), we find that a(r) presents 
an initial contracting phase {h(0) < 0) provided that Wo < 1. Since viable models imply 
an expanding Universe at late times (what is guaranteed by Eq. |26] once W increases and 
becomes greater than unity), models with Wo < 1 necessarily have a minimum in the scale 
factor and, hence, they are non-singular. 

In the opposite, when Wo > 1, the Universe evolution starts with an expanding phase 
(h > 0) which will continue at any later time. As a matter of fact, the requirement x(r) > 
together with Eq. ( |14D imply that /(r) > y(r)/2 for any r > 0. Consequently, h{r) is 
always positive. 
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2. Early-time solutions 



Using /o — y /2 the first-order approximation of Eqs. ( fiSp and (|49|) reduce to 

u = y (55) 
| [7i - l2/a]r (a > 72/71) 



(56) 



^-vr(^) (a = 72/71 



where 



a > 72/71 (57) 

is required according to Eq. fl35]). 

Consequently, when a > 72/71, the early-time behavior of the dynamical fields and 
variables is 

2 2 

$ oc r^i Q -T2 oc t 3 "- 1 (58a) 

q — 1 1— a 

a oc ttl q -t2 oc ti^ta (58b) 

! 71Q ,_(2-3 7 ) 37 1 _ n 

j 2[ 7ia - 72 ] OC t 2 l-3a ( a jL 1) 

(58c) 
const. > (a = 1) 

l-3a 

t oc r ti"-t2 (58d) 

while, when a = 72/71, we find that solutions have exponential behaviors which are not 
compatible with a coupling function 3 + 2co>($) as that assumed in Eq. (j38"l) . 

Since a must satisfy Eq. (|57|) , the above solutions imply that $(0) = 0. Consequently, 
by identifying a 2 = (k + 1)/A 2 and e = 2/(71 a — 72), the coupling function (p9|) can be 
expressed in the form given by Eq. fl38"D . 

The scale and speed-up factors have instead a more complicated behavior depending on 
the a value: 

i) If a > 1 (k;(0) > 0), the scale factor vanishes at r = while £ diverges to +00. Models 
are then singular. 



17 



ii) If a < 1 (ll>(0) < 0), then a(0) diverges to +00 and £ diverges to —00. Models are 
then non-singular. 

Hi) a — 1 (^(0) = 0), the scale and speed-up factors have a finite non- vanishing value at 
t = while $(0) = 0. In this case, the a and $ derivatives are both positive at r = and, 
hence, a negative r value exists where the scale factor has a singular point and $ < 0. We 
then note that our choice of the origin of times ensures that gravitation is always attractive 
for t > 0. However, in this case, an alternative choice of the origin of times could also be 
the requirement a(0) = (implying a repulsive gravitation at early times). 

All these behaviors are shown in Figures 2, which represent the analytical solutions 
obtained when the coupling function is given by (3 + 2u)/3 = a 2 + br (a = 1) (see appendix 
B), as well as the numerical solutions obtained for (3 + 2u)/3 = a 2 + br 1 ! 2 (a ^ 1) . We 
also see from these figures that models with a < 1 can present a non-monotonic behavior of 
the speed-up factor. 

C. Models with and y = ($'(0) > 0) 

1. General features at arbitrary times 

According to Eq. (25), the expanding or contracting behavior of the scale factor, a(r), 
is determined by the sign of [/ — yj (2y/W)]. 

Taking into account that W(t) is a monotonic increasing function (l/yW < 1/y/Wo) 
and that integration constants must satisfy fo = yo/2 = 0, Eq. (|16D implies that y{r) < 
7 2 t/VW / o. Consequently, y/VW < 72T/W0 an d, using Eq. (|TTD (with f = 0), we obtain 

[/ - y/(2VW)} > (71 - 7 2 /^ )r/2 (59) 

which implies two different situations in the time evolution of a(r), depending on the initial 
value of W(t): 

i) Wq > 72/71: in this case, [/ — y/(2y/W)] > for any r > 0, and there exists a 
undefined expansion. 
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ii) Wo < 72/71: A similar argument as that used in the paragraph preceding Eq. (|59"D 
implies that for any time previous to that, r 1; in which W reaches the value W\ = 72/71: 

[/ - y/(2VW)] < (60) 

Consequently, there exists an initial phase of contraction in the scale factor (for r < r{). 
Since f(r) > y(r)/2 (x(r) > 0), Eq. (25) implies an expansion process when W becomes 
larger than unity. Then, the scale factor must present a minimum (Fig. 3b), in the time 
interval between W = 72/71 and W = 1. 



2. Early-time solutions 

According to Eqs. fl3~4|) and Q35|) , integration constants must satisfy f = y /2 = and 
a > 72/71- Consequently, the u and t> functions can be approximated by 

u = ^(71 + 7 2 /a)r (61) 



I (71 - l2/a)r (a > 72/71) 
(« = 7 2 /7i) 

When a > 72/71, the early-time behavior of models is then given by: 



(62) 



$ oc r 4 ^ /A oc t 472/B (63a) 

[ r~ 2 ^ r *m/fl ^ 72 / 7l ) 
i oc ^ (63c) 

[roci (a 2 = 72/71) 

t oc r B / A (63d) 

A = 7l V - 7 2 2 ; 5 = 2 72 + 3 7 ( 7 ia 2 - 72 ) (63e) 

while, when a = 72/71, solutions are exponentials which are not compatible with Eq. (|38|) 
We then find that all models imply $(0) — > 0. Consequently, by identifying a 2 = 
(k + 1)/ X 2 and e = 472/ (7 2 a 2 — 7I), the coupling function ( |39"D can be expressed in the form 
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given by Eq. (p8|) . The scale and speed-up factors have instead a wider variety of early- 
time behaviors. If a 2 > 72/71 (e < 2), models are singular (a(0) — > and £(0) — > +00). 
If a 2 < 72/71, models are non-singular with a(0) — > +00 and £(0) — » —00. Finally, if 
a 2 = 72/71 (e = 2), models are again non-singular but now with a(0) ^const. (a minimum) 
and £(0) — > 0. These behaviors are shown in Figures 3, which represent the analytical exact 
solutions for W = a 2 + br (case a 2 > 72/71 of appendix B), and the numerical solutions for 
W = a 2 + br 2 (a 2 = 72/71) and W = (a + qr 3 ) 2 (a 2 < 72/71). 

D. Models with a = 

1. General features at arbitrary times 

According to Eq. (34) these models are only compatible with the initial condition /o = 
— 2/0 > 0. Consequently, $(r) is a monotonic decreasing function and Eq. (p6|) implies that 
the scale factor a(r) increases indefinitely from a vanishing value (x(0) = 0) at the origin of 
times. This behavior is shown in Figures 4, which display the numerical solutions obtained 
when (3 + 2u)/3 = 6r(l + cr 1/2 ) 3 (b, c> 0) at any time. 

2. Early-time solutions 

When 3 + 2uj vanishes as r — > 0, the coupling function can be expressed as 

1 3 + 2uj I /3 = 6r e (64) 
with 6 > and < e < 2. In this case, the first-order approximation of u and v functions 



will be given by (see Eqs. [5(J and |T|) : 



u = r( 2 ^)/ 2 (65a) 

>/6(2 - e) 

u = = -Vo (65b) 

Then, Eqs. @-(H) give 
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$ OC T 



2/(72+5) 



t -2/(l+5) 



(66a) 



a oc r 



1/(72 +<5) 



OC t 



1/(1+5) 



(66b) 



72+2 (j-l) 72 +2(g-l) 
£ OC T 2 (72+'5) OC t 2( 1 + 5 ) 



(66c) 



t oc t^+s 



(66d) 



where 



72e 



(67) 



2 - e 



corresponds to the exponent of |$ — 1| appearing in Eq. (p8|). 

As expected, we find that all models are singular and $(0) — > +oo. The early behaviour 
of £ depends instead of the 5 value. If 5 > 5 C = (37 — 2)/2, we find that ^(0) — > +00. If 
5 = S c , ^(0) has a finite positive value. Finally, if 5 < 5 C , £ vanishes at r = and presents a 
non-monotonic behavior (see Figures 4). In any case, since $ diverges to infinity at r — >• 0, 
we can approximate $ ~ $ — 1 and Eq. ( plj ) becomes the coupling function given by Eq. 



According to Eq. (|i~4"D, models with W < are characterized by the fact that x(t) = 
P + 2/ 2 /4- Consequently, 2 has always a positive value which can be initially vanishing only 
if integration constants are chosen so that they satisfy the condition given by Eq. (p6|): 
fo = Vo = 

Since the factorization P0| ) is not further possible, solutions to the scalar-tensor cos- 
mological equations can only be obtained by applying the general procedure described in 
Eqs. (^2])-(0). The early-time expressions can be deduced by considering Eq. (|39| ) and 
approximating y{r) by the first non- vanishing term of: 



(§|) with k = 0. 



V. MODELS WITH 3 + 2lu < 




(68) 
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where integration constants must be chosen according to the criteria established in Section 



m~B 

A. Models with a ^ and y < 0, ($'(0) < ) 
1. General features at arbitrary times 

Although we already know that all these models have a stationary point at r = (see Eq. 
p7|) , we will now analyze the conditions needed so that such a stationary point corresponds 
to a minimum in a(r). 

From the stationary point condition (Eq. f = y / (2yJ\W \)} < 0, the function /(r) 
is initially negative, but it becomes positive for any r > t\ = 2|/ |/7i- Consequently, we 
can ensure that the scale factor expands for r > T\ (see Eq. ^) and, therefore, it presents 
a minimum value at r = provided that: 

/' = 7 i/2 > W{2y/\W\)Y (69) 

If the early-time behavior of W(t) is represented by Eq. (|39|), the above condition is 
equivalent to e > 1, with the additional constraint b\yo\/2 < a^ia 2 + 72) for e = 1. 



2. Early-time solutions 

Since x(0) = fo + Uo/^ > 0, condition ( pTj) implies that a(0) and $(0) are finite and 
non- vanishing, while £(0) = 0. Consequently, as a first-order approximation: 

$ — > const. > 1, a — > const. > 0, £ — > 0, , t oc r (70) 

Using Eq. (^) to obtain a higher approximation, we find: 

$ ~ $(0)(1 + cir) (71a) 

a~a(0)(l + c 2 r 2 ) (71b) 

£ oc r (71c) 

t oc r (71d) 
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where a(0) > 0, $(0) > 1 and 

4a 

ci = n - ^ 2 (72a) 

_ 7ia 2 + 72 + y b /2a 

c 2 = 2 , ~ 2 ^ , (e = 1) 72b 

2/5(1 + a 2 ) 

Models are then non-singular and the scalar field takes a finite non-vanishing value at 
early-times (see Figure 5). It must be noted that, from Eq. ( |71a|) , the only e value consistent 
with Eq. (BSI) is e = 1. 



B. Models with a ^ and y > 0, ($'(0) > 0) 

1. General features at arbitrary times 



According to the conditions established in Sect. [HFE, $(t) is a monotonic increasing 
function with an initial value < $(0) < 1 . 

On the other hand, since /(r) and |W(r)| are monotonic increasing functions while y(r) 
decreases monotonically, Eq. ( pTf ) implies 

[/ - y/(2^iwi)] > [/o - y /{2j\Wo\)] = (73) 

and, according to Eq. (|26|) , a{r) then expands at any time from an initial minimum value 
o(0) > (see Fig. 6) 

2. Early-time solutions 



By following the same procedure as in Sect. [V A| , we obtain again the solutions ([H]), but 
now with < $(0) < 1. 

C. Models with a ^ and y = 0, ($'(0) < ) 
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1. General features at arbitrary times 



From the initial condition (^) , / = y = 0, we find that /(t) > and y{r) < for 
any r > (strictly). Consequently, $(r) decreases monotonically while, from Eq. (p6|) . 
a(r) expands indefinitely This behavior is shown in Figure 7, which displays the numerical 
solutions obtained in models defined by |W(r)| = a 2 + br at any time. 



2. Early-time solutions 
Using /o = yo = 0, the y(r) function can be approximated (see Eq. |68[) by 



y(r) = -(72/«)r (74) 



Thus, we obtain 



$ oc r- 4 ^/ A oc t- 472/B (75a) 

a oc r 2 (7i« 2 +72)/A K t 2( 7l a2 +72)/B (75b) 

e oc r 272/A oc t 272/B (75c) 

t oc r fi / A (75d) 

A = 7l 2 a 2 + 7 2 , 5 = 3 77l (l + a 2 )-8 (75e) 

Models are then singular (with a = a' = 0) and $(0) — > +oo. 

It must be noted that, for 7 values implying B < 0, the above early-time solutions 
correspond to the limit t — > —00 instead of t — > 0. This is also true for any other asymptotic 
solution where i is proportional to a negative power of r. 



D. Models with a = 

7. General features at arbitrary times 



As we have seen in Sect. D ITB, when 3 + 2u;(0) vanishes from negative values, the only 
possible choice of integration constants is f = y = ($'(0) — 0)- J us ^ like in the previous 
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subsection, this choice ensures that /(r) > and y(r) < for any r > 0. Consequently, 
$(t) is a monotonic decreasing function while, from Eq. (p6|), o(r) expands indefinitely from 
an initial singularity (x(0) = 0). This behavior is shown in Figures 8 through the analytical 
exact solutions obtained when |W(t)| = br e at any time (see also Appendix B). 



2. Early-time solutions 



Using Eq. (|68|) with fo = Uo = 0, the y(r) function can be approximated by: 

272 _(2- e )/2 



» (T) = -<^wi T (e < 2) (76) 



Thus, we obtain: 



$ oc r" 4 / A oc t" 4/B (77a) 

a oc r 2 ' A oc t 2 / B (77b) 

i oc r 3 * 1 -^ oc t 2 ^- 5 )/ 5 (77c) 

t oc r B / A (77d) 

A = 7 2 + 25, fi = 2(5-7 2 + 2 (77e) 

where 5 = 726/ [2 (2 — e)] is the exponent appearing in Eq. (|38~1) . 

Models are then singular and the scalar field decreases from $(0) — > +00. The initial 
value of the speed-up factor depends instead on the 5 value. If 8 > 1 , £(0) diverges to +00, 
(Fig. 8a); if 8 = 1, £(0) has a finite non- vanishing value; finally, if 5 < 1 , £(0) becomes zero, 
(Fig 8b). 

Table 1 summarizes all the early-time asymptotic solutions that we have found in sections 
[TV}-|V| for 7 = (false- vacuum models) and 7 = 1 (matter-dominated models). This table 
also contains the early-time solutions for vacuum and radiation-dominated models (7 = 4/3) 
obtained in our previous work ||15| . The way in which these results can be applied to derive 
the asymptotic solutions of other scalar- tensor theories (not strictly defined by Eq. |38[) is 
illustrated in the Appendix A. 
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VI. LATE TIME SOLUTIONS 



We will analyze in this section the way in which scalar-tensor theories converge towards 
GR at late times. We note however that, although solar-system experiments require in fact 
that any viable scalar-tensor theory must converge towards GR during the matter-dominated 
era, such a condition is not compulsory for earlier epochs in the universe evolution. Serna Sz 
Alimi [T21J have in fact shown that some scalar-tensor theories are able to predict the right 
primordial abundances of light elements even when such theories are very different from 
GR during primordial nucleosynthesis. Consequently, at the end of the inflationary epoch, 
a scalar-tensor theory can deviate considerably from GR and still be compatible with all 
astronomical observations. 

When convergence to GR (3 + 2u — ► ±00) is imposed at late-times (r — > 00), Eqs. (|12|), 



([14]) and fll7|) imply y' — > and x(t) ~ jfr /4. The integration method described in Sect. 



Ill then lead to solutions which, as a first-order approximation, are similar to those found 



in GR: 

$ ~ f ~ 1 (78a) 

a oc r 2/71 (78b) 

I ln(r) (7 = 0) 
t oc I (78c) 

|^ r 37/7i ( 7 ^ 0) 

A higher-order approximation to the $ and £ solutions can be found by considering (see 
Eq. [H) that, in the limit $ — ► 1, the asymptotic form of the coupling function is 

\3 + 2lu(<!>)\/3 = \- 2 /\<!>-l\ 5 = br e (79) 

with A > 0, 5 > 0, b > 0, and e > 0. 

Using Eq. (|7jJ), the y{r) function is given by 

y(r)=y(r 1 ) + -^- ) (r-r 1 )^; (e ^ 2) (80) 
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where T\ denotes the time at which Eq. (|79|) starts to be a good approximation for the 
coupling function. The e = 2 value in Eq. fl80| ) has been excluded because it leads to models 
which are incompatible with Eq. fl38|). 

In the limit r — > oo, we can find two possibilities for the y function: a) y(r) converges 
towards a finite value, y c , and b) y(r) diverges to ±00. We will now analyze separately the 
solutions corresponding to these two cases. 



A. Case y — ► y c , (1 < 5 < 2) 



Since $ ~ 1, y{r) ~ y c and x(r) ~ 7 2 r 2 /4, Eqs. ([TO]) , ([XT]) and (|?J) imply: 



(81) 



which, together with Eq. (P7|), leads to: 

$ = 1 - sign(y c )cr" 2/(2 " 5) (82a) 

e = l + ab c 2 /(2 7l 2 )]r- 2 (82b) 

where 1 < S < 2 (e > 2) and 

5 = 2e/(2 + e) (83a) 

A" 2 = 6c 5 (83b) 

c ^ [2A(2 - 5)|2/ c |/ Tl 2 ] 2/(2 - 5) (83c) 

We then find that <3> converges towards the GR value from below ($ — > 1~) or above 
($ — ► 1+) depending on the sign (positive or negative, respectively) of y c . On the other 
hand, £ converges to unity from above (£ — > 1 + ) or below (£ — > 1~) depending on the sign 
(positive or negative, respectively) of 3 + 2u. 



B. Case y -> ±00 (6=1) 

In this case, $ ~ 1, y(r) oc 2cr 72 /[&(2 - e)]r (2_e)/2 (see Eq. gD[) and x(r) ~ 7 2 r 2 /4. 
Consequently, Eqs. (fLOf), (|Tl"l) and (|79|) imply: 
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$ = 1 - acT" e (84a) 
i = 1 - <x(c/4)T- £ [(e - l)/(2 - e)](4 - lx e)r^ (84b) 

with 5 = 1, 1/X 2 = be, and 

6 = 1 ± [1 - (8A 2 72 /7i)] 1/2 ; (8A 2 < 71V72) (85) 

We then find that $ converges towards the GR value from below ($ — > 1~) or above 
($ — > 1+) depending on the sign (positive or negative, respectively) of 3 + 2u. On the other 
hand, £ converges to unity from below (£ — > 1~) or above (£ — > 1 + ) depending on the sign 
(positive or negative, respectively) of <r(e — 1)(4 — 7ie). 



VII. CONCLUSIONS 

We have presented in this paper a method to derive exact solutions for K = scalar- 
tensor cosmologies with an arbitrary u;($) function and satisfying the general perfect fluid 
state equation P = (7 — l)pc 2 , where 7 is a constant and < 7 < 2. This procedure 
has some common aspects with the indirect method previously proposed by Barrow and 
Mimoso jnj]. In particular, solutions are not directly found through the coupling function 



cu($) which defines each specific scalar-tensor theory, but by means of a 'generating function' 
and a suitable change of variables. A given choice of the generating function produces, after 
completely solving the field equations, a particular form of co , ($). Therefore, the specific 
class of scalar-tensor theories under study is not known 'a priori', but only 'a posteriori'. 

Unlike the previous indirect method, the non-linear transformation of variables used in 
our 'semi-indirect' procedure allows us to take the time dependence of the coupling function 
itself as the generating function. This has the obvious advantage that much of the main 
properties characterizing the functional form of uj(Q) can be easily known 'a priori', what 
results specially useful to analyze the asymptotic behavior of scalar-tensor cosmologies. 

Using this procedure, we have supplied a comprehensive study of asymptotic cosmological 
solutions in the framework of a wide class of homogeneous and isotropic theories. We 

28 



have also described the qualitative behavior of models at intermediary times and, for some 
particular scalar-tensor theories, we have obtained solutions which are exact at any time 
characterized by a constant 7 value. Our analysis then covers all the main epochs in the 
cosmic history (inflationary, radiation- and matter-dominated models) and, therefore, it 
extends other works (specially those of Refs. |TJj and |T7|] ) also devoted to a systematic study 
of scalar-tensor models. All the early-time solutions found in this paper are summarized in 
Table 1. We find that singular models with 3 + 2uo > and yo 7^ have always early-time 
solutions for a and $ which are independent of 7. Consequently, in these models, all the 
epochs in the universe evolution have the same early behavior. 

We finally note that the solutions obtained in the present paper, together with those 
of other previous works [|15|,[17]], enable complete cosmological histories to be constructed 
through initial vacuum (or false vacuum), radiation, matter, and final vacuum-dominated 
eras. The wide diversity of possible scalar-tensor cosmological models can be restricted from 
additional constraints like those obtained from the primordial abundance of light elements 
Z0|| , but also from their implications on the initial spectrum of density perturbations []2T 



and other strong-field tests (e.g., Ref. |2"2" 
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APPENDIX A: CONNECTION WITH OTHER SCALAR-TENSOR THEORIES 

The asymptotic solutions found in this paper are not restricted to scalar-tensor theories 
with a coupling function 3 + 2uj strictly given by Eq. (^) . Such asymptotic solutions hold 
for any arbitrary form of 3 + 2uj admitting (at early or late times) a series approximation as 
that expressed by Eq. fl3"E|). 

Let us consider, for example, a scalar-tensor theory defined at any time by 
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with B > and $ e (0, 1). At early times ($ — > 0), we can write 1/(1 — $) ~ 1 + $ and, 
hence, the Taylor approximation of Eq. (|86|) is 

3 + 2^ = ^[1x1(2) - 1+ 1 1 - $ r 1 ] (87) 
In (2) 

2/ 



which has the form given by Eq. (|38| ) with A 2 = 3 In (2)/B, k = ln(2) — 1 and 5=1. 



Consequently, 1/a = \/^/k + 1 = y31n(2)/£>, and we see from Table 1 or Eqs. ( [58a| )-( |58"rj|) 
that its early solutions at any epoch are 



ct-l 2 



aocts^-i; $ a t 3 «- 1 (88) 

Models are then singular provided that a > 1 (B > 31n(2)) while, otherwise, they are 
non-singular. 

At late times ($ -> 1, 7 = 1), we have In [($ + l)/2] ~ -(1 - $)/2 and 

3 + 2^ = 25 (89) 

which is again of the form given by Eq. (||) with A 2 = 3/(25), k = and 5 = 1. According 
to Eqs. (p4|)-(p5|), this theory converges at t — > 00 towards the GR solutions provided that 
A 2 < 9/8, that is, provided that B > 4/3. 

In the same way, much other theories (as the late-time behavior of all the classes of 
theories defined in [O) are asymptotically described by the solutions obtained in this paper. 
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APPENDIX B: ANALYTICAL EXACT SOLUTIONS 



The exact solutions obtained by applying the method described in this paper to those 
theories defined by (3 + 2w)/3 = a 2 + br e (a > 0, b > 0) are: 
1) If f Q = y /2 > 0, e = 1, v = 0, and a = +1 



$ OC 



z — a 



z + a - 272/71 



— 4 atan f 214+52' 



£71^ + 72 



v vtrt y i? < 



i? = 



1+2/y/fl 1+2/-/H 

(2 — zi) ti (2 — z 2 ) 71 -R > 



a-l a-(2-3 7 )/ 71 

a oc (2 — a)Ti«-T2 [2; + a — 272/71] Ti a -T2 < 



where 



M i/7i e 7 1 ^ a tan(2i£+22 ) J R<0 (90b) 



yl/7lg -71(71^+72) 



i2 = 



(90a) 



u = (7i/2)r + (72/6) [Va 2 + 6r - a] + w 

2 = V« 2 + br 
-72 + 

Zi = 



7i 



-72 

^2 = 



v^R 



7i 



i? = [71a + 72] 2 - 2kjiU 



2) If y Q = 0, e < 2, a 2 = and a 

, 2(2-6) 

$ = [1 + Cr^]~ 



a = ( 7l r/2) 2 / 7l [l + CY- £ ] ™ 

7ie-j 'Yi e — 4 

^72 



7ie- 4 



T- e [l + CT- e }~ 



where 



C = 



72 



(91a) 
(91b) 
(91c) 

(91d) 
(91e) 

(92a) 
(92b) 
(92c) 

(93) 



b V 7i(2-e)y 

It is straightforward to see that, in the limits r — > and r — > 00, the above expression 
reduce to the corresponding asymptotic expressions of Sections III-V. 
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TABLE 1 asymptotic solutions at early times: \3 + 2u\^ (3/\ 2 )(k+ | $ - 1 \~ 5 ) 
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Vacuum and Radiation-dominated models 
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FIGURES 




FIG. 1. Models with W > 0, a / and j/o < 0: numerical solutions obtained by considering 
the theory defined, at any time, by W(t) = (p + qr) 3 with b = 1, 7 = 1, and a) a 2 = 0.2, yo = —3, 
b) a 2 = a 2 = 1/9, 7/0 = — 3, c) a 2 = 0.01, yo = — 1- Solid lines represent the scale factor, dotted 



lines represent the scalar field, and dashed lines represent the speed-up factor. 
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FIG. 2. Models with W > 0, a / and j/o > 0: matter-dominated (7 = 1) solutions obtained 
by considering the theory defined, at any time, by a) W{t) = a 2 + br 1 ^ 2 with a 2 = 25/9, 6 = 1, 
y = 1, b) W(t) = a 2 + br with a 2 = a 2 c = 1, 6 = 1, y = 1, c) a 2 = 4/9, 6 = 0.5, y = 0.02. 
Symbols for lines are the same as in Figure 1 
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FIG. 3. Models with a / and yo = 0: matter-dominated solutions obtained by considering 
the theory defined, at any time, by b = 1 and a) W(t) = a 2 + br (a 2 = 5/9), b) W(t) = a 2 + br 2 
(a 2 = a 2 = 1/3), c) W(t) = (q + qr 3 ) 2 (b = 2aq and a 2 = 7/27). Symbols for lines are the same 
as in Figure 1 
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FIG. 4. Models with a = and IF > 0: numerical solutions obtained by considering the theory 
defined, at any time, by W(t) = 6r(l + cr 1 / 2 ) 3 (b, c> 0) with y = -2 and a) 7 = 1 (5 > 5 C = 1/2), 
b = 0.1, and c = 10, b) 7 = 1.12 (<5 < 5 C = 0.68, 6 = 0.1, and c = 10, c) 7 = 10/9 (5 = <5 C = 2/3), 
6=1, and c = 0.5. Symbols for lines are the same as in Figure 1 
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FIG. 5. Models with W < 0, a ^ and yo < 0: numerical solutions obtained by considering 
the theory defined, at any time, by W(t) = a 2 + br with 7 = 1, a 2 = 1, b = 1, yo = —2. Symbols 
for lines are the same as in Figure 1 
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FIG. 6. Models with W < 0, a / and yo > 0: numerical solutions obtained by considering 
the theory defined, at any time, by W(t) = (p + qr) 3 with 7 = 1, a 2 = 1, b = 10, yo = 1- Symbols 
for lines are the same as in Figure 1 
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FIG. 7. Models with W < 0, a ^ and yo = 0: numerical solutions obtained by considering 
the theory defined, at any time, by |W(r)| = a 2 + br with 7 = 1, a 2 = 1, b = 1. Symbols for lines 
are the same as in Figure 1 
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FIG. 8. Models with W < and a = 0: analytical solutions obtained by considering the 
theory defined, at any time, by |W(r)| = br € with b = 1, 7 = 1 (so that the critical value S c = 1 
corresponds to e c = 4/3) and a) e = 1.2, b) e = 1.4. Symbols for lines are the same as in Figure 1 
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